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Abstract: This research work appertains to the acceptance sampling plan under the neutrosophic 
statistical interval method (ASP-NSIM) based on gamma distribution (GD), Burr type XII 
distribution (BXIID) and the Birnbaum-Saunders distribution (BSD). The plan parameters will be 
determined using the neutrosophic non-linear optimization problem. We will provide numerous 
tables for the three distributions using various values of shape parameters and degree of 
indeterminacy. The efficiency of the proposed ASP-NSIM will be discussed over the existing 
sampling plan in terms of sample size. The application of the proposed ASP-NSIM will be given 
with the aid of industrial data. 
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1. Introduction 


In most industries, an acceptance sampling plan (ASP) plays a vital role in testing the quality of 
the manufactured product. The ASP in statistical quality control has established its importance to 
decide on the acceptance or rejection of a lot of a product. At the stage of inspection, it may not 
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possible to inspect a lot of the product. In such a case, a sampling plan is a very appropriate 
alternative for deciding whether a lot of the product should be accepted or rejected based on 
verification of the quality in a selected sample. Sampling plans are essential when lot sizes are very 
large, the chance of assessment errors is high, and the cost of the inspection is extremely significant. 
These ASPs were first applied at the time of the Second World War by the Americans. [1] developed 
an ASP and promoted its application in the industry. The authors said that “it is a method of 
inspecting a sample of products to decide whether the product lot is to be accepted or not, based on 
the results obtained and if the number of failures during the test time does not exceed the acceptance 
number then the lot is accepted.” Hence, the ASP is the best indispensable method for the 
checkup/testing of the product. The plan parameters which are used in the checkup/testing of the 
manufactured goods are ascertained in accordance with the specified producer’s risk as well as the 
consumer’s risk. For more details, see [2]. Thus, the well-defined sampling plan will minimize the 
risks and provide a smaller sample size for testing the manufactured goods. More details on the ASPs 
from truncated life tests could be encountered in [3] for log-logistic distribution, [4] for Weibull 
distribution for repetitive group sampling, reference [5] reliability sampling plans based on hybrid 
censoring, [6] for reliability acceptance sampling plans for the exponential distribution, reference [7] 
for Birnbaum-Saunders distribution. 

The aforementioned ASPs are classical sampling plans which are developed for different 
statistical distributions under the postulation of the known quantities. All of the ASPs mentioned 
above are only useful when the researcher confirms the percentage of nonconforming items in the 
manufactured goods. In recent years, the fuzzy approach has become a popular research 
methodology in the field of ASPs for uncertainty in the fraction of nonconformity. Various 
researchers boosted their research contribution on the design of ASPs based on the fuzzy atmosphere 
for different situations comprising as [8] studied a design for a single sampling attribute plan based 
on fuzzy sets theory, [9] developed the sampling plans by attributes using the fuzzy approach, [10] 
presented the fuzzy acceptance sampling plan and its characteristic curves, [11] studied single 
sampling plan with fuzzy parameter, [12] developed a single acceptance sampling plan with a fuzzy 
parameter for Poisson distribution, [13] proposed the double sampling plan using fuzzy Poisson 
distribution, [14] developed the sequential sampling plan using fuzzy SPRT, [15] proposed 
inspection error and its effects on single sampling plans with fuzzy parameters, [16] further studied 
the acceptance sampling for the influence of TRH using crisp and fuzzy gamma distribution and [17] 
studied a fuzzy mathematical analysis for the effect of TRH using acceptance sampling plans. 

Because of its versatility in addressing the uncertainty in the data and the parameters of the 
distributions, neutrosophic statistics (NS) or neutrosophic logic has drawn the attention of more 
researchers during the past several years. [18] developed neutrosophic statistics and which is the 
generalization of traditional statistics which can be employed to analyze the data in uncertain 
situations. [19] studied a new attribute ASP-NSIM and pointed out that “in the case of uncertainty, 
the existing sampling plan cannot be applied for the inspection of a lot of the product due to in a 
practical situation it is not necessary that under some circumstances all the observations/parameters 
are determined values”. For more information about the NS and its applications, please refer to [18]. 
The analysis of neutrosophic numbers from rock measurements is studied by [20] and [21]. In recent 
times, [22] developed different sampling plans using NS. Reference [23] suggested using a 
determinate sample size for the sampling plans. 

To the best of our knowledge, there is no work based on ASP-NSIM for the GD, BXHD and 
BSD. In the present article, we will develop an ASP for the GD, BXIID and BSD under the NISM. 
For the proposed ASP-NSIM the neutrosophic plan parameters will be determined for various 
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degrees of indeterminacy. Extensive tables will be presented for practical use. The application of the 
proposed ASP-NSIM will be given with help of industrial data. 


2. Some continuous distributions under a neutrosophic environment 
2.1. Gamma distribution 


Let Y,, €[Y,,¥, ];i=1,2....,m, denote the neutrosophic random variable where Y, refers to the 
lower value of the indeterminacy interval and Y,, presents the upper value of indeterminacy interval. 
The neutrosophic form of Yyi is Yui = Yii + Yuilyy; lyy€lly,» ly yl. where Y,; is the lower value of 
a random variable and Yyjly, is the indeterminate (upper) value of the random variable and 
Iyy€[ly,» ly, is the measure of the degree of indeterminacy. 

The neutrosophic cumulative distribution function (NCDF) of GD is given below: 


1 
BAY G00 =F b, 2 Vy 20,5, € [BeBe oy e lost| (1) 
Td, o 


N 


where tom is the lower incomplete gamma function, b, e[b,,b, | is the neutrosophic 
N 


shape parameter and oO, € epee | is the neutrosophic scale parameter. The neutrosophic form of 
the shape parameter is by = b;, + bylyy; Loy €Uby Hoy | and the neutrosophic form of the scale 
parameter is Oy = 0, + dylan; eles ja where the first value shows the determinate value of 
shape and scale parameters and the second value shows the indeterminate parts and I, yell re ly, | 
and Toy € loys Ioy| are the measure of degrees of indeterminacy, respectively. 

For the present study, it is assumed that the neutrosophic shape parameter b, is known. When 
by is unknown, it can be estimated from the available data. The average lifetime of neutrosophic 
GD is fy =byo, . A product failure probability before the time Jy) is denoted as 
Py =F (Yy <Yyo) and is defined as 


r= pe r[du2 (2) 


Here, we express neutrosophic termination time ),y. is a product of constant a and neutrosophic 
mean life fly, 1.€., Vo =ALMyy- Therefore, Eq (2) could be rewritten in terms of neutrosophic mean 


Hy as follows: 


i allno 
= ——¥(b ) 
by 7 Ly /by 
1 
— Y(by,aby / 2°). (3) 


N 
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2.2. Burr-type XII distribution 


In the area of quality and reliability, the Burr-type XII distribution is the more intriguing distribution. 
Additional examples of this distribution's use in reliability analysis and quality control may be found in [24]. 
Consider a lifetime of the product follows a neutrosophic BXIID with NCDF is given below: 


‘N 


ky | 
Poni) =1-[ 1422) | Vn 20, A, > 0, ky > 0,V x, >0 
(4) 
whereas, 2, €[A,,2,] is the neutrosophic scale parameter, ky €[k,,k,]and vy €[v,,v, | are 


neutrosophic shape parameters of neutrosophic BXIID. The neutrosophic form of the shape 
parameter is ky = ky, + kylky; Ney € icy» Ley | and the neutrosophic form of the scale parameter is 
Vy = Vy + Vylyy; Ing €[Ty,) lal where the first value shows the determinate value of shape and scale 
parameters and the second value shows the indeterminate parts and Trey € [Ley I ku | and Iyy€ [Lop joe 
are the measure of degrees, respectively. 


The median lifetime of the product using neutrosophic BXIID is givenas wy =A, (2" "y — 1)" ; 


A product failure probability before the time V,,. is defined as 


_ Sih by YN 
Py =1 (2) | : (5) 


It is convenient to determine the neutrosophic termination time ),) as a multiple of the 
specified median life Jd), that is Vy =ALy . Hence, the neutrosophic failure probability given in 


Eq (5) can be expressed as follows: 


a(2" Si N 
Dy =1—| 1+ (6) 
Hy | Hyo 
2.3. Birnbaum-Saunders distribution 
The NCDF of the neutrosophic BSD is given by 
1 % 5 A 
F(yy) =®y | — (22) (%) Vy 2 UO ee 0, SO, (7) 
N one Yn 


where 7, €[7,,7%,] is the neutrosophic shape parameter, 5, €[5,,6,] is the neutrosophic scale 
parameter of NBS distribution and ®,,(.) is NCDF of standard normal distribution. The neutrosophic 


El, where the first value shows the 


Ti: 


form of the shape parameter is yy = Y, + Vulyyilyy€[y 


determinate value of shape and the second value shows the indeterminate parts and I, yell 
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2 
The mean of neutrosophic BSD is given by mw, = dy [i+22). A product failure probability 


before the time ,, is denoted by p, and it is given below: 


hs 
1 
Py =®,|— [2 [2] : (8) 
Vy \\ Oy no 
The neutrosophic termination time Vy. is expressed aS Vy) =Aflyy. The neutrosophic failure 


probability given in Eq (8) could be expressed as 


_ sacl wn 
i a,[ 1428) / 
Hy | H 
Py =®,| —4| ————- | - aa ; (9) 
Yn Hy | Ho t ) 
Ay| i+ 5 


3. Methodology of acceptance sampling plans under NSIM 


The proposed ASP-NSIM is stated as follows 
Step1. Select a random sample of size n from a yielded lot and conduct the life test for sample 
items for the specified time yy. 
Step2. Count the number of items that failed before the termination time y,, and denoted as d. 
Step3. The lot could be accepted when d is less than or equal to c before termination time yy 
otherwise reject the lot. 


We are interested to determine the plan parameters (n,c) of the proposed ASP-NSIM through 
neutrosophic nonlinear optimization for which given, producer's risk (@ ) and consumer's risk ( £ ) 


are satisfied. 

The plan parameters (n,c) will be determined using a two-point approach, as the producer 
always wishes that the probability of acceptance of the lot should be greater than 1-q@ at an 
acceptable reliability level (ARL), say p,, and the consumer wishes that the probability of 
acceptance should be smaller than f at the lot tolerance reliability level (LTRL), say p,,. The 


neutrosophic operating characteristic (NOC) function will be used to determine the plan parameters. 
The NOC of the proposed ASP-NSIM is given below: 


Po(0n)= 3" abt (10 


The plan parameters of the proposed plan can be determined through the following neutrosophic 
non-linear optimization 
Minimize nn. (11) 


Subject to 
e(n 3 
Po (0s)= (4 pilD) *21-a (12) 
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Py (P02) = Johalt— Dal <B n2=2,c20. (13) 
d=0 

To find the plan parameters, the neutrosophic non-linear optimization is implemented as follows: 

Step 1. Fix the values of a, 6, shape parameter of specified distribution, a, degree of indeterminacy. 

Step 2. Determine the values of n and c such that the above-mentioned constraints are met. 

Step 3. There are many combinations of n and c that satisfied the given constraints. Choose that 

combination of n and c where n is minimum. 

To obtain the plan parameters of the proposed ASP-NSIM, the quality of the product is 
constituted in terms of the ratio of true mean/median life and specified mean/median life. In this 
study, the quality level is measured through the ratio of its mean lifetime to the true mean lifetime, 
Ty = Hy / Myo - To find the plan parameters of the proposed plan, the values ofa and # are fixed in 
advance. As mentioned earlier, the producer requires the lot acceptance probability should be larger 
than 1—a@ at ARL, say p,,. The producer is interested to have sampling plan parameters at 


various values Ny = Ly / Myo = 2,4,6,8,10,12, on the other hand, the consumer is interested to have 
the plan parameters at My / Ay) =1. 


The plan parameters plan parameters (n,c) are determined using the above-mentioned 
neutrosophic non-linear optimization. The plan parameters (n,c) for the three distributions are 
reported in Tables 1-3. The plan parameters n and c are chosen such that n is minimum for a 
given degree of uncertainty. The plan parameters n and c for GD when a@ =0.05; £ =0.25, 0.10, 


0.05; % =2,4,6,8,10,12; a=0.9, 1.0, 1.1 and by=0.10 are given in Table 1. The plan parameters n 

and c for BXIID when a@ =0.05; £ =0.25, 0.10, 0.05, 0.01; 7“ =1.2,1.5,1.6,1.8,2,4; a=0.9, 1.0, 1.1; 
ky=2.0, Vy=0.5 are given in Table 2. The plan parameters and c for BSD when a@ =0.05; ££ =0.25, 
0.10, 0.05, 0.01; % =2,4,6,8,10,12; a =0.9, 1.0, 1.1 and yy=3 are given in Table 3. 


From Tables 1—3, we noticed the succeeding points: 
(1) For fixed @ and shape parameters for various distributions the values of n and c 
decreases as ratio 7%, increases. 


(2) For other same parameters, the values of mand c decrease when the # increases. 
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=0.10. 


Table 1. Neutrosophic plan parameters for GD when by 
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Table 2. Neutrosophic plan parameters for BXI[D when ky=2.0, vy=0.5 and Ixy = Ivy = In- 


a=0.90 a=1. 0 a=1.10 
B | ty | ty =0 | Iy=0.05 | Iy=0.10 | Iy=0.20 | Iy =0 | Iy=0.05 | Iy=0.10 | Iy=0.20 | Iy =0 | Iy =0.05 | Iy =0.10 | Iy =0.20 


n fe {|n {ce f{[n fe f{n fe {n fe f{n fe {[n fe {n fe {n fe fn c n c n c 
1.2 | 90 | 33 | 76 | 27 | 69 | 24 | 55 | 18 | 91 | 37 | 78 | 31 | 69 | 27 | 53 | 20 | 93 | 41 | 80 | 35 | 67 | 29 | 57 | 24 
1.5 | 33 | 10 | 28 26} 7 | 21/5 | 30] 10] 28/9 | 23 |7 | 21 |6 | 30] 11 | 28 | 10 | 23 | 8 21 | 7 
1.6 |20|5 | 18 18 | 4 15 | 3 18 | 5 | 18] 5 16 | 4 13 | 3 19 |6 | 17 |5 14 | 4 12 | 3 

Mee 1.8 | 15 | 3 12 | 2 12 | 2 13 | 3 | 13 | 3 13 | 3 11 | 2 15 12 | 3 12 | 3 10 | 2 
2.0 | 12 12|}2 |9 1 11 11 | 2 11 |}2 |8 1 12|3 | 10 | 2 10 | 2 7 1 
4.0 | 12 | 2 1 1 2 1 1 10 1 1 
1.2 | 67 | 25 63 
1.5 | 23 

ART 1.6 | 16 


10 
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a=0.90 a=1.0 a=1.1 
B Ny | ly Tyy 9.05 | Ty, =0.10 | ty, =0.20 | Ly, =0 Iyy =9.05 | 1y,=0.10 | Iy,y=0.20 | 1, =0 Tyy =9-05 | Ty =0.10 | ly, =0.20 
n c n c n c n c n c n c n c n c n c n c n c n c 
2 246 | 138 | 234 | 130 | 224 | 123 | 196 | 105 | 246 | 144 | 231 | 134 | 221 | 127 | 199 | 112 | 242 | 147 | 229 | 138 | 221 | 132 | 199 | 117 
4 88 | 46 | 80 | 41 | 77 | 39 | 72 | 35 | 86 | 47 | 80 | 43 | 77 | 41 | 72 | 37 | 83 | 47 | 82 | 46 | 79 | 44 | 70 | 38 
6 51 | 25 | 46 | 22 | 45 | 21 | 40 | 18 | 49 | 25 | 49 | 25 | 46 | 23 | 42 | 20 | 52 | 28 | 49 | 26 | 46 | 24 | 42 | 21 
ae 8 35 | 16 | 35 | 16 | 32 | 14 | 29 | 12 | 35 | 17 | 34 | 16 | 30 | 14 | 27 | 12 | 37 | 19 | 32 | 16 | 33 | 16 | 30 | 14 
10 | 28 | 12 | 26 | 11 | 24 | 10 | 23 | 9 26 | 12 | 27 | 12 | 23 | 10 | 24 | 10 | 25 | 12 | 27 | 13 | 24 | 11 | 23 =| 10 
12 | 22 |9 20 7 
189 | 107 | 182 | 102 | 173 | 96 | 155 | 84 | 188 | 111 | 176 | 103 | 169 | 98 | 153 | 87 | 186 | 114 | 176 | 107 | 171 | 103 | 155 | 92 
68 | 36 | 65 | 34 | 62 | 32 | 54 | 27 | 65 | 36 | 64 | 35 | 61 | 33 | 55 | 29 | 66 | 38 | 63 | 36 | 62 | 35 | 56 | 31 
6 38 | 19 | 37 | 18 | 37 | 18 | 32 | 15 | 40 | 21 | 35 | 18 | 37 | 19 | 33 | 16 | 40 | 22 | 37 | 20 | 34 | 18 | 31 | 16 
a 8 28 | 13 | 24 | 11 | 26 | 12 | 21 |9 28 | 14 | 27 | 13 | 23 | 1 | 22 | 10 | 27 | 14 | 27 | 14 | 24 | 12 | 23 | Il 
10 | 22 | 10 | 21 |9 19 | 8 17 |7 21 | 10 | 22 | 10 | 20 |9 18 | 8 22 | 11 |21 |10 |19 |9 18 | 8 
12 | 19 | 8 17 |7 15 | 6 16 | 6 18 | 8 16 | 7 14 | 6 15 | 6 17 | 8 17 | 8 17 | 8 14 | 6 
114 | 66 | 108 | 62 | 104} 59 | 90 | 50 | 111 | 67 | 107 | 64 | 103 | 61 | 91 | 53 | 112 | 70 | 108 | 67 | 99 | 61 | 89 | 54 
40 | 22 | 37 19 
aioe 6 23 | 12 | 25 ll 
8 16 | 8 18 7 
10 | 13 | 6 11 6 
12 | 9 4 9 4 
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4. Comparison study 


In this section, we will study the efficiency of the proposed ASP-NSIM for the GD, BXIID, and 
BSD distributions with ASPs under classical statistics. The proposed ASP-NSIM will be more 
flexible and informative than the traditional sampling plans, see [21]. For a fair comparison, we will 
consider the same values of specified parameters for the proposed ASP-NSIM and the existing 
sampling plans under classical statistics. The values of n for the proposed ASP-NSIM and the 
existing sampling under classical statistics for GD are shown in Table 4. The values of n for the 
proposed ASP-NSIM and the existing sampling under classical statistics for BXIID are shown in 
Table 5. The values of n for the proposed ASP-NSIM and the existing sampling under classical 
statistics for BSD are shown in Table 6. A plan which provides smaller values of the sample size is 
known as an efficient sampling plan. From Tables 4-6, it can be noted that the proposed ASP-NSIM 
provides smaller values of n as compared to the sampling plans under classical statistics. Moreover, 
the existing sampling plan does not yield information regarding the measure of indeterminacy. Thus, 
we conclude that the proposed ASP-NSIM GD, BXIID and BSD are more efficient than the sampling 
plan under classical statistics and capable of handling the uncertain situation. 


Table 4. Comparison of sample sizes in the proposed plan and the existing plan for GD 
when /,,,=0.20. 


— a=0.90 can 
B Bayo Proposed Classical Proposed Classical 
5 0 55 22 23 
4 9 22 10 Me 
0.25 | 6 6 13 2 : 
g 6 ll 5 5 
104 9 : 2 
35 45 33 38 
4 16 18 14 A 
0.10 | 6 Ul 1 9 a 
8 9 p 8 
10 6 9 5 : 
5 0 2%6 22 23 
4 9 1 10 me 
0.05 | 6 6 9 3 : 
8 6 6 5 5 
10 4 6 : 2 
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Table 5. Comparison of sample sizes in the proposed plan and the existing plan for 
BXIID when Iy = 0.20. 


a=0.90 a=1.0 
P neue BrOpONe Classical plan Proposed Plan Classical plan 
Plan 
1.2 55 90 53 91 
LS 21 33 21 30 
0.25 | 1.6 15 20 13 18 
1.8 12 15 11 13 
2 9 12 8 11 
1.2 42 67 39 70 
15 16 23 17 26 
0.10 | 1.6 11 16 9 14 
1.8 8 10 9 12 
2 8 10 i 9 
1.2 28 4] 27 42 
15 11 15 10 16 
0.05 | 1.6 8 10 q 10 
1.8 6 5 7 
2 6 S) 5 


Table 6. Comparison of sample sizes in the proposed plan and the existing plan for BSD 


when Lyy = 0.20. 


a=0.90 a=1.0 
B Bane Proposed Plan Classical plan pes Classical plan 
Plan 
2 196 246 199 ae 
4 nN 88 72 ae 
0.25 |6 40 51 - - 
8 29 35 27 — 
1023 28 24 a8 
2 155 189 153 ioe 
4 54 68 55 oe 
0.10 32 38 2S - 
1 28 22 28 
10 17 22 18 zal 
90 114 91 111 
3] AO 31 40 
0.05 21 23 a a 
13 16 14 ms 
10 13 13 Ul 2 
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5. Exemplification of the proposed plan 


In this section, the application of the proposed ASP-NSIM will be discussed. [25] discussed that 
the lifetime is imprecise in practice. [26] used the lifetime data in intervals. More information on 
imprecise lifetime data can be seen in [27]. Assume that a producer wants to supply the mean life 
assertion for his manufactured goods, and he asserts that the true mean life of the manufactured 
goods is /4y)= 1000 hours. The quality supervisor judges after verification of the product and decide 
whether the producer’s claim for the lifetime of the manufactured goods is valid or not. Assume that 
the test time of the experiment is y,,)= 1000 hours. Thus, we get the termination ratio for the 
experiment is a=1.0. Considering the consumer’s risk is as #=0.25, the producer’s risk a =0.05, 
the ratio of its mean lifetime to the true mean lifetime, dy / Hy  =2, and the shape parameter of the 
GD is taken as by=0.10. By following [25], we assume a degree of uncertainty that is /,,=0.20. 
Using this information, from Table 1, the values of the plan parameters of the proposed ASP-NSIM 
are n=22 and c=12. It shows that the required sample size for testing a lot of the product is 22. 
Suppose the quality supervisor has picked out 22 samples from a lot of a product and the test is 
conducted up to the specified time of 1000 hours. The submitted lot of the product is accepted if no 
more than 12 items failed before the time 1000 hours. The proposed ASP-NSIM can be applied to 
other distributions on the same lines. 


6. Conclusions 


In this article, we developed the ASP-NSIM for the GD, BXIID and BSD. The plan parameters 
for the proposed ASP-NSIM were determined using neutrosophic non-linear optimization. Extensive 
tables were presented for practical use at various degrees of uncertainty. From the comparative study, 
it can be concluded that the proposed ASP-NSIM is more efficient than the sampling plans under 
classical statistics. The proposed test using other statistical distributions can be extended as future 
research. The proposed plan using a cost model can be studied in future research. 
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